Formation of matter-wave soliton molecules 
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We propose a method of forming matter-wave soliton molecules that is inspired by the recent 
experiment of Dris et al. [lj. In the proposed set-up we show that if two solitons are initially 
prepared in phase and with a sufficiently small separation and relative velocity, a bound pair will 
always form. This is verified by direct numerical simulation of the Gross-Pitaevskii equation and 
by the derivation of the exact interaction energy of two solitons, which takes the form of a Morse 
potential. This interaction potential depends not only on the separation but also on the relative 
phase of the solitons and is essential for an analytical treatment of a host of other problems, such 
as the soliton gas and the Toda lattice of solitons. 
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Introduction. - Interactions between solitons have 
been the subject of intense study in the past few decades 
@, & y< This was motivated both by the fundamen- 
tal physics involved and by the prospects of important 
technological applications [5j. At the fundamental level, 
solitons show a dual nature by behaving as classical par- 
ticles but also exhibiting their wave nature when they, 
for instance, scatter off each other. In particular, it was 
first pointed out by Karpman and Solov'ev using per- 
turbation theory Q , by Gordon using the exact two soli- 
tons solution Q , and by Anderson and Lisak employing a 
variational approach [4| , that at large distances the force 
between two solitons decays exponentially with their sep- 
aration and is proportional to the cosine of their phase 
difference. From the applications point of view, this in- 
teraction between solitons is actually a problem since it 
leads, when using optical solitons in fibers as data car- 
riers, to the destruction of the information stored in a 
sequences of solitons. 

The exciting possibility of forming so-called soliton 
molecules gives yet another stimulus to the problem. 
Such molecules have indeed recently been realized ex- 
perimentally with optical solitons by Stratmann et al., 
who also measured the relative phase of the two solitons 
and studied their binding mechanism Matter- wave 
soliton molecules have, however, not yet been created ex- 
perimentally. Theoretically, soliton molecules have been 
studied previously in two-dimensional Bose-Einstein con- 
densates without and with Q dipolar interactions. In 
the latter case an effective molecular soliton-soliton po- 
tential was found, which results from the dipolar nature 
of the condensate and interlayer effects. Experimentally, 
bright solitons and soliton trains were observed as the 
remnants of a collapsing condensate that was created by 
using a Feshbach resonance to switch the interaction from 
repulsive to attractive Q. It was found that neighbor- 
ing solitons repel each other in spite of the fact that the 
interatomic interactions arc attractive. A variational cal- 



culation showed that the repulsive force is caused by a 
7r phase difference between the solitons [lfj , in agree- 
ment with known result for optical solitons [2, 5, 4]. The 
conditions in these experiments are thus not favorable 
for forming bound states of solitons. Instead, if the two 
solitons start with a zero phase difference, they attract 
each other and may ultimately form a molecule. Our 
analysis of the soliton-soliton interaction energy shows 
that while this is indeed a necessary condition, it is not 
sufficient for soliton molecule formation. The initial sep- 
aration and relative velocity need to be sufficiently small 
such that the interaction is not dramatically weakened by 
the exponential tail and the kinetic energy of the relative 
motion does not considerably exceed the soliton-soliton 
interaction energy. 

An adjusted version of a recent experiment by Dris 
et al. [l| is suggested here as a possible experimental 
method of realizing soliton molecules in strongly elon- 
gated Bose-Einstcin condensates. In the experiment of 
Ref. [H, a trapped bright soliton was launched onto a 
potential barrier at the center of a harmonic trap. Here 
partial transmission and reflection created two solitons, 
which then after half an oscillation in the trap recom- 
bine again into a single one. We will show, however, that 
a soliton molecule can be formed if both the harmonic 
trapping potential and the potential barrier are switched 
off at the time when the two solitons reach their classical 
turning points. This guarantees the relative velocity to 
be zero. Moreover, the two solitons are also guaranteed 
to have a zero relative phase then, since they were cre- 
ated coherently from the same soliton. The separation 
between the solitons prepared in this way is controlled 
by the curvature of the harmonic trap. Thus, it appears 
that all the conditions for soliton molecule formation can 
be met in such an experimental set-up. In the following, 
we show that this is indeed the case via direct numer- 
ical simulations of the Gross-Pitaevskii equation. This 
is then followed by an analytical analysis of the soliton- 
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soliton interaction energy, which also confirms the fea- 
sibility of our proposal. Moreover, finding an accurate 
interaction potential for two solitons at an arbitrary sep- 
aration, is important not only for the present problem of 
forming soliton molecules, but may also lay the grounds 
for a host of other problems involving solitons interac- 
tions such as the soliton Toda lattice and the soliton 
gas [Hj]. 

Soliton molecule formation. — The evolution of soli- 
tons in the original experiment of Ref. as well as in 
our modified proposal, can be described by the effectively 
one-dimensional Gross-Pitaevskii equation (GPE) 



9(z,t)=0, (1) 



where, here and throughout, lengths are scaled to the 
characteristic length a z — ^JK/rnuj z of the harmonic po- 
tential mcu^z 2 /2, time to 1/uj z , and the wavefunction 
ty(z,t) to 1/ y^2a z u!±/ui Zl with cu z and tu± the axial and 
radial frequencies of the strongly elongated (u± 3> oj z ) 
trapping potential, respectively. In these units, the 
strength of the interatomic interaction will be given by 
the ratio g = —a/a z , where a is the s-wave scattering 
length. A square or gaussian potential barrier is added 
to the harmonic potential at the center that is included 
in Vext(z). 

We now show that Eq. ((T|) indeed accounts for the in- 
teresting recombination of the two returning solitons at 
the center, as can be seen in Fig. HJa). Repeating this 
simulation for the case that we switch off the harmonic 
trapping potential and the potential barrier at the point 
when the two solitons reach their classical turning points, 
results in the formation of a molecule of two solitons, as 
shown clearly in Fig. [ljb). The relative intensities or 
densities of the two solitons in the molecule are set by 
the initial velocity of the single original soliton and the 
height of the potential barrier. 

In another version of this experiment, one can also 
think of phase imprinting one of the two solitons once 
it reaches the classical turning point. In such an experi- 
ment, the phase of the cycle performed by the molecule is 
controlled. For instance, one can imprint a 7r phase dif- 
ference leading initially to an increase of the separation 
of the solitons. If their initial separation is small enough, 
however, they will ultimately return back and form the 
molecule since their relative phase will change from n to 
zero as they move away from each other and the force 
becomes attractive, as will be shown next. 

Soliton- soliton interaction energy. — As mentioned in 
the introduction, the interaction energy of two solitons 
depends on their separation and on their relative phase. 
It is also known that such an interaction may lead to the 
formation of a bound state of two solitons. The key point 
here is that dynamically the relative phase of two solitons 
changes as their separation changes. If two solitons, that 
are initially in phase, approach each other sufficiently 
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FIG. 1: (Color online) Spacio-temporal density plots for (a) 
the experiment of Ref. [l| and (b) our proposal for soliton 
molecule formation. In (a) a single soliton is split by the 
potential barrier into two solitons that upon returning to 
the trap center recombine into the original one. In (b) the 
harmonic trapping potential and the potential barrier, shown 
with the solid blue curve, are switched off at t — 15.75 when 
the two solitons reach the classical turning point. Initial soli- 
ton parameters: velocity = —0.8, amplitude = 0.8, position 
= 5, and interaction strength g = 0.0113. Square barrier's 
parameters: height = 0.75 and width = 0.5. 



slowly, their relative phase changes such that the force 
between them becomes repulsive. When the two solitons 
then move away from each other, the force becomes at- 
tractive again and the cycle completes itself. Thus, the 
manner by which the relative phase changes with the sep- 
aration between the solitons is essential for understand- 
ing their binding properties. This fact is qualitatively 
understood in the literature, but has never before quanti- 
tatively been captured in a formula that accurately gives 
the interaction potential of the solitons for both short as 
well as large separations. Here, we arrive at this objec- 
tive by deriving the exact dependence of the phase, and 
hence the interaction potential, on the solitons separa- 
tion. In addition to being of importance by itself, this 
potential provides a convenient analytical understanding 
for the possibility of soliton molecule formation. 

Using the Inverse Scattering Method, the exact two 
soliton solution of the homogeneous version of Eq. (JXJ) is 
derived in Ref. In Ref. a re-derivation puts the 
solution in terms of 8 parameters, namely the two initial 
center-of-mass positions Zx,2 of the solitons, their center- 
of-mass velocities fi,2, their phases $1,2, and their square 
amplitudes ni,2- The separation between the solitons can 
then be extracted analytically, from which also the force 
follows in terms of these parameters. Here, we extend 
this formalism to derive the interaction energy and to 
show that it takes the form of a Morse potential. 

For simplicity, the relative velocity is set to zero and 
the limit of nearly equal soliton densities n_ = ni — n\ <C 
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n + = ri2 + ri\ is taken. The separation A(t) between the 
solitons and their relative phase 3>(t) turn out to be given 
by 



A(t) = - log 
a 



a (a 2 + 2a(3 cos (w i + $ ) + /3 2 ) 
(3 g 2 n 2 _ 



(2) 



and 



d>(t) 



tan 



tan" 



/3sin (w i + 0o ) 
a + /3 cos (wt + 0o ) 
1 / a sin i + 0o ) 
/3 + a cos + 0o ) 



(3) 



where fj = exp (a xo/4), oj = agn^/8, and a = gn + . 
In addition, xq and 0o are two constants that can be 
determined from the initial conditions A(0) = Ao and 
<fr(0) = <f>o- The equations of motion for A and $ 
can now be obtained by differentiating twice with re- 
spect to t. The resulting equations are simplified by 
inverting Eqs. (|2|) and (J3j> to express cos(wi + o ) and 
sin (ut + 0o ) in terms of A and Finally, we obtain in 
this manner 



A 



a 
32 



COS<I>, 



sin $, 



(4) 



(5) 



which are identical to Gordon's formulae |3j. Here, how- 
ever, we have the exact solutions of these equations of 
motion, namely Eqs. ([2} and (j3|). This allows for ob- 
taining the force A in terms of A by solving Eq. ([2j for 
cos (u) t + 0o ) and substituting the solution in Eq. (QJ. In- 
tegrating the resulting expression with respect to A leads 
to the soliton-soliton potential energy V ss = — J A dA 



V s , 



(cu 2 /4) cos$ 



/ e -ia(A-A. q 



2e" 



i Q (A-A oq ) 



(6) 

which has the form of a Morse potential with an equilib- 
rium position at A eq given by 

A A 4 i ^ 2 -2a/?cos$o + /? 2 ^ 

Here, (3 — /3(Ao, $0) corresponds to the real solution of 
the initial conditions. 

Alternatively, the potential can be determined from 
the exact solution directly. To that end, the location of 
the two solitons are calculated numerically from the den- 
sity \^(z, t)\ 2 , from which the solitons separation and the 
associated acceleration are computed. Integrating the re- 
sulting force numerically, we again obtain an interaction 
potential. In Fig. [2j we plot the interaction potential 
obtained by the two methods. It is clear from this fig- 
ure that the numerics captures the main features of the 



potential in Eq. © apart from a slight change in the 
location of the equilibrium point and depth of the po- 
tential. The change is due to the fact that the difference 
in the center-of-mass positions is numerically difficult to 
define uniquely and is therefore not exactly equal to the 
parameter A in the exact solution when the two solitons 
are too close together and interference effects take place. 
The shape of this potential, specifically the existence of 
a global minimum, proves that a bound state of two soli- 
tons can be formed provided the initial separation and 
relative phase are such that the energy is negative. 

Finally, we notice that both the equilibrium point and 
depth of the potential depend on the initial conditions. 
For further clarification of this rather unusual situation, 
we plot in Fig. [3] the solitons separation A(t) and rela- 
tive phase $(t) during one period of the molecule's os- 
cillation. It is clear that by increasing the initial sepa- 
ration, the equilibrium point around which the solitons 
oscillates also increases. Notice that, for both cases, the 
phase changes from to n over half the period 7r/w of 
the molecule's oscillation, which means that the force 
changes from attractive to repulsive over the same pe- 
riod of time. On the other hand, the dynamics of A(t) 
depends on Ao since for larger Ao there is initially a 
smaller overlap between the solitons and thus they will 
approach each other more slowly. By the time the force 
has changed sign, the solitons of the molecule with larger 
Ao have a final separation that is larger than the one 
with smaller Ao , but both have now to start moving away 
from each other, which amounts to different equilibrium 
points. We have also checked that the same conclusions 
are drawn directly from the exact solution of the GPE. 




FIG. 2: Interaction potential between two bright solitons in 
terms of their separation. Solid curve is obtained numerically 
from the exact two solitons solution of the Gross-Pitaevskii 
equation. Dashed curve corresponds to formula ((6]). Parame- 
ters: n_ = 0.2, n+ = 4.7, g = 0.5, $0 = 0. 

Other applications. — Having an analytic formula for 
the interaction potential as a Morse potential, it is ap- 
pealing to considering further applications. In particular, 
we can consider a classical gas of solitons interacting with 
the Morse potential in Eq. (|6]) . Here, the thermodynam- 
ics is expected to be enriched by the dependence of the 
potential on the relative phase of the solitons. 

The equation of state of such an imperfect gas 
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FIG. 3: (Color online) Two upper red curves: solitons sep- 
aration A(t). Two middle green curves: phase difference 
&(t). Two lower blue curves: force A(t). Time range cor- 
responds to one period of the molecule's oscillation. Two 
values of the initial separation were used: Ao = 24.0 (solid 
curves) and Ao = 24.5 (dashed curves). Parameters: $o = 0, 
n_ =0.1, n+ = 1, g = 1. 

can be derived by a virial expansion P/ks T = p + 
£?2 p 2 + S3 p 2 + ■ ■ ■ , where p is the number density 
and B2 and B3 are the second and third virial coeffi- 
cients. The second virial coefficient is given in terms 
of the two-body partition function which reduces to 
the following integration on the soliton-soliton potential: 
-(1/2) J [exp {-V ss /k B T) - 1] dA. Due to the fact that 
V ss depends on the initial conditions and in the partition 
function all configurations should be taken into account, 
an additional integral on all possible initial conditions 
should be performed. A simpler but equivalent approach 
is using the potential that corresponds to Eq. (j4|, namely 
V ss = — (ct 2 /2) exp (—a A/4) cos<E>, where $ is treated as 
an independent thermodynamic variable and the initial 
conditions are absent from the potential. In the limit 
of weak interactions a -C fcsT, a simple expression can 
be derived B 2 = -(ir/8)a 3 / (k B T) 2 , which leads to the 
equation of state 

P 7T a 3 2 

Exploring the possibility of a phase transition to a soliton 
liquid requires extending the virial expansion to the third 
virial coefficient, which involves three-body processes and 
hence requires the knowledge of the three solitons solu- 
tion, which is beyond the scope of this Letter. 

As another interesting application, we briefly consider 
a Toda lattice of solitons, which was successfully used to 
model the dynamics of a train of A^-solitons in external 
potentials [14J. The force on the i-th soliton at position 
Zi will be given by 

z i = F ss (A i+1 )-F S s(A i - 1 ), (9) 

where A i+ i = z i+1 - Zi, = %i - Zi-i, and F ss = 

—dVss/dA. The ground state of such a system will be a 
series of equidistant solitons with peak intensities alter- 
nating between two slightly different values. The value of 



n_ cannot be arbitrarily small since the nearest-neighbor 
separation between the solitons in the lattice diverges 
as logn_ [l3[|. On the other hand, n_ should be small 
enough that the solitons are sufficiently separated and 
that the coalescence condition n + g = 1 is avoided [l3| ■ 
In this configuration the net force on any soliton will 
vanish and the lattice moves as a rigid object that can 
be considered as a "soliton of solitons" . Collective exci- 
tations and sound waves can be considered as the low- 
energy excitations. 

In conclusion, we have shown that a molecule of two 
bright solitons can be realized in attractive Bose-Einstein 
condensates. This was demonstrated both by numerical 
simulations and by showing that the interaction poten- 
tial between the two solitons is of a molecular type and 
takes the form of a Morse potential. Applications of this 
potential in the problems of soliton gas and Toda lattice 
of solitons were pointed out. 
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